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PACS 75 . 47 . Lx - Manganites 

pacs 75 . 47 . Gk - Colossal magnetoresistance 

PACS 72 . 80 . Ng - Disordered solids 

Abstract. - Disorder on the active d element site is usually very disruptive for conduction 
and long range order in perovskite transition metal oxides. However, in the background of phase 
competition such 'B site' dopants also act to promote one ordered phase at the expense of another. 
This occurs either through valence change of the transition metal or via creation of 'defects' in 
the parent magnetic state. We provide a framework for understanding the complex variety of 
phenomena observed in B site doped manganites and identify the key parameters that control the 
physics. Using a spatially resolved analysis of B ions in various manganite phases we explain the 
existing data and predict new situations where highly polarisable phase separated states can be 
created. 
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Correlated electron systems like the cuprates and man- 
ganites involve competition between various long range or- 
dered phases [1,2]. The interplay of this phase competition 
with weak disorder underlies phenomena like cluster coex- 
istence, percolative transport, and colossal response. The 
nature of disorder seems to be crucial for these effects, as 
observed in the manganites [2], and a low concentration 
of impurities on the active d element site is an effective 
trigger for phase separation [3-12] and the associated per- 
colative effects. 

The results of such substitution depend on the reference 
state and the chemical nature of the impurity. For man- 
ganites, with rare earth (RE) and alkaline earth (AE) com- 
bination REi-a^AEa^MnOs, several intriguing results exist 
for Mn site doping, (i) Magnetic dopants like Cr [3-7] , Co 
or Ni [8] (but not Fe) on the Mn site in a x = 0.5 charge or- 
dered insulating (CO-I) manganite promote a percolative 
ferromagnetic metal (FM-M) , while non magnetic dopants 
of the same valence do not. {ii) The orbital ordered A type 
antiferromagnet (AF) at x = is destabilised in favour of 
a ferromagnetic state [9] by both magnetic and non mag- 
netic dopants, {in) In contrast to the cases above, where 
charge-orbital order is suppressed^ doping Fe on a ferro- 
magnetic metal [10, 11] at x ~ 0.4 promotes a charge or- 
dered insulating state! On spatial imaging most of these 
systems reveal phase separation (PS) and many of them 
also exhibit enormous magnetoresistance. It is vital to 
uncover the organising principle behind this diversity of 



effects, if we are to exploit B site disorder as a tool for 
phase control. 

There is unfortunately no microscopic model, let alone a 
theory, for randomly located B dopants in the manganites. 
In this paper we write down the first detailed model for 
B impurities in a manganite host, and study the effect of 
these dopants in a variety of manganite phases using a real 
space Monte Carlo technique. 

Our principal results are the following: (i) We discover 
that the following hierarchy of effects arise in all B doping 
cases: (a) change of the effective valence on the Mn sites, 
(b) percolation of the metallic phase through impurity free 
regions, and (c) 'reconstruction' of the background mag- 
netism and charge order by magnetic dopants, (ii) By 
exploring the prominent manganite states, and different 
B dopants, we are able to explain most of the outstanding 
experimental results, (iii) We suggest a new experiment 
to test out an unexplored insulator- metal transition driven 
by B site disorder, (iv) We demonstrate how B impurity 
locations determine the percolation pattern and may allow 
atomic level control of current paths in a material. 

The simplest classification of B site dopants is in terms 
of their valence in the manganite host. Among the usual 
dopants Zn, Mg, and Co, are divalent, z.e, in a 2+ state, 
Ni, Cr, Fe, Sc, and Al are trivalent, while Ru, Sn, and 
Ti are tetravalent. Some elements can exist in multiple 
valence states, e.g^ Ni can also be +2, and Ru can be +5, 
but that will not affect our qualitative arguments. The 
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valence, a, of the dopant affects the effective carrier den- 
sity on the Mn sites through the charge neutrahty require- 
ment on the compound RE^t^AE^+Mn^l^^B^Og", where 
T] is the % of B site doping, and we write the Mn valence 
as 3 + i^. This yields iy{r]^ a, x) = {x -\- 77(8 — a)) /(I — 77). 
The effective Cg electron count on Mn is n = 1 — mod- 
ified from no = 1 — xat7^ = 0. This change of effective 
carrier density can itself drive phase change as we will see 
later. Secondly, dopants with same valence can have dif- 
ferent effects depending on their magnetic character. Non 
magnetic dopants only affect the Mn valence, while those 
with partially filled d shells can have magnetic coupling to 
the neighbouring Mn moments. Experiments suggest that 
Cr has strong AF coupling [13] to the Mn ions, Ni couples 
ferromagnetically [13], while Fe, despite its magnetic 
configuration, couples rather weakly. 

Based on the inputs above we construct the following 
model for manganites with B site dopants: Hfot = Href + 
Himp H~ Hqoup') with 

{ij)a i {ij) 

i i 

Himp = y ^liaa^naa 
naa 

a[3 

Hcoup = ^ ^a^/S^liaa^jPa + S^.Sj + Vc QnQj 

{nj)a {nj) {nj) 

The reference 'manganite model' Href is constructed to 
reproduce the correct sequence of phases in the 'clean' 
limit. It involves the nearest neighbour hopping of Cg 
electrons with amplitude [14], Hund's coupling Jh, 
AF superexchange J between Mn spins, and Jahn- Teller 
(JT) interaction. A, between the electrons and the phonon 
modes Q^. The stiffness of the JT modes is K. We will 
not consider RE-AE cation disorder in Href since such 'A 
site' disorder masks the B site effects. The sites Ri,Rj 
and operator c^c'^ refer to Mn locations, and S^, etc, are 
Mn spins. The local physics of the B ions is contained in 
Himp where refers to the B locations and the operators 

d^ refer to the B ion Cg states at an energy V above the 
center of the Mn band. The sites R^ are random, with 
only the constraint that two B dopants are not nearest 
neighbours, to minimise electrostatic repulsion. Hcoup in- 
volves {i) Cg hopping matrix elements t^fd^ which we keep 
the same as between the Mn, (ii) for magnetic B ions, a 
superexchange coupling J' between the B moment and 
the neighbouring Mn moments, and {in) a nearest neigh- 
bour Coulomb repulsion Vc between the B dopant and the 
neighbouring Mn. The total charge qj on the Mn ion is 
4 — nj , where nj is the Cg occupancy, and Qn is the (fixed) 
B ion valence. Fig.l is a schematic, showing the rele- 
vant levels on Mn and B, and the coupling between these 
atoms. Earlier attempts at modelling B dopants (or Mn 
vacancies) only employed a random onsite potential [15]. 
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Fig. 1: Colour online: Energy levels of the B and Mn ions and 
the couplings between Mn-Mn and Mn-B. We show a schematic 
density of states at the left highlighting the primarily Mn band 
and the broadened B level. 

The overall carrier density is controlled through 
the chemical potential term: -K^iaa 4aa^iaa + 
^naa^nacrdnacr)' We usc the Standard limit Jn/t > 1, 
and set K = 1. In studying magnetic field effects we will 
use a coupling Hmag = -^-{^i ^i^^n ^n), where h = zh 
is the applied field. We treat all spin and phonon degrees 
of freedom as classical [16], and measure all energies in 
units of the Mn-Mn hopping t [14]. The spins, S^, etc, are 
treated as unit vectors, and the magnitude of the spin is 
absorbed in the couplings J and J^ 

The parameter space of the problem is very large. We 
have to contend with manganite states at different doping 
(x), different (inverse) bandwidth (A/t), and AF strength 
{J/t). The impurities are characterised by their concen- 
tration (t^), the Cg level (V/t), and the magnetic cou- 
pling (J'/t). We use X/t = 1.6 as a typical JT cou- 
pling, J/t = 0.1 (to capture the correct phase competition 
around x ~ 0.5), and Vc = O.lt [17]. We will explore sev- 
eral doping levels x = 0.25, x = 0.4 and x — 0.5 and study 
the effect of a low density of B dopants. 

We use a Monte Carlo (MC) technique based on the 
'travelling cluster approximation' (TCA) [18]. It allows 
ready access to system size ~ 40 x 40 using a moving 
cluster of size 8x8, and handles disorder accurately. 
The method, and the associated transport calculation, has 
been extensively benchmarked by us and used successfully 
in several earlier studies [14,19]. 

The first principle that we wish to illustrate is phase 
separation driven by change in effective carrier density. 
Fig. 2 shows the phases in Href foi" varying 'hole density' 
X — 1 — n. The T = phases, in increasing order of x, are 
an orbital ordered insulator at x = [19], a ferromagnetic 
charge ordered insulator (FM-CO-I) at x = 0.25, a FM-M 
window between x ~ 0.30 — 0.42, the CE-CO insulator at 
X = 0.50, and a magnetic phase ('A-2D'), with structure 
factor peaks at q = {0, tt}, {tt, 0}, between x ~ 0.55—0.60. 
Between these phases are the shaded windows of PS. If the 
carrier density is in one of these PS windows the system 
would break up into coexisting patches of the two adjacent 
phases. For a system at the edge of PS a small valence 
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Fig. 2: Colour online: The x — T phase diagram of our refer- 
ence model in 2D at X/t = 1.6 and J/t = 0.1. The true 3D 
transition temperatures can be estimated roughly as 3/2 times 
the 2D scales indicated here. The phases of interest to us are 
the ferromagnetic CO insulator at x = 0.25, the ferromagnetic 
metal for x - 0.3 - 0.42, and the CE-CO phase at x = 0.50. 
Shaded regions indicate phase separation. 



change driven variation in the carrier density can push it 
into the PS window. Fig. 2 shows many such possibilities. 
These can be exploited by choosing dopants of suitable 
valence. The tendency towards large scale PS competes 
with the fragmenting effect of disorder, leading finally to 
a percolative state. 

Let us start with x = 0.25. In the clean limit the orbital 
ordered (00) JT insulator at x = [19] is separated from 
the FM-CO-I at X = 0.25 by a wide PS window. The 
FM-CO-I can be looked upon as the charge ordering of 
doped holes with double the lattice periodicity in both x 
and y directions. In Fig. 2 the phase is bounded on the 
right also by a PS window, which separates it from the 
homogeneous FM-M. To 'metallise' the x = 0.25 phase we 
need impurities that decrease the effective electron density, 
pushing it towards the FM-M. From our expression for n = 
1 — u at low T] we obtain n ~ (1 — ^) —r]{3-\-x — a)-\-0{r]'^)^ 
so any impurity with valence < 3 will serve to lower n 
and push the system into the x ~ 0.25 — 0.30 PS window. 
The system could then phase separate, with FM-CO-I and 
FM-M clusters whose pattern is controlled by the B ion 
locations. 

Fig. 3 top row shows a snapshot of the charge density 
field rir in the reference CO state at 77 = (left), followed 
by the impurity locations (center), and Ur (right) in the 
presence of the B ions. The impurities are non magnetic 
(J' = 0), divalent, with V = 5 (which is roughly half the 
Mn bandwidth) and r] = 2%. Trivalent impurities yield 
similar results. The emergence of homogeneous regions of 
n ~ 0.7 {x ~ 0.3) coexisting with short range (SR) CO 
patches is clear. The FM-M live in the 'impurity free' re- 
gions, while the SR-CO can coexist with the impurities at 
short distance but loses coherence at larger scales. The 
percolation of the FM-M patches creates a 'global' metal- 




Fig. 3: Colour online: The density field rir at low temperature. 
Columnwise: left - Ur in the reference state (77 = 0), center - 
location of the B ions, right - Ur in the presence of the B ions. 
The top row is for the x = 0.25 CO-I state, with 2% doping 
of a non magnetic 2+ ion with 1/ = 5. Middle row is for the 
FM-M at X = 0.40 with 6% doping of a non magnetic 2+ ion 
with y = 5. Bottom: doping on the CE-CO-I at x = 0.50, 
with a magnetic 3+ ion: 77 = 6%, 1/ = 5, = 0.2. 

lie state. Fig. 4. (a) show the evolution of the resistivity 
p{T) with T]. Since both the reference CO phase and the 
emergent metal are FM we cannot get any magnetoresis- 
tance in this regime. It would be interesting to choose y 
in the (Lai_2yPr2y)o.75Cao.25Mn03 family so that the ma- 
terial is at the T = insulator-metal phase boundary, and 
explore the impact of 2+ or 3+ B dopants. 

Let us shift to x = 0.4 which is next to the PS window 
separating the FM-M from the x = 0.5 CE-CO-I phase. 
The reference state is a homogeneous metal. We explore 
the impact of non magnetic divalent dopants with V = 5 
on this state, and push n into the PS window. Such a sit- 
uation had been explored experimentally early on by dop- 
ing Mg [12] onto PrcyCa^ySrca-^yMnOa. A metal-insulator 
transition in the ground state was seen with increasing Mg 
doping. Similar results were observed on doping Fe [10] 
(which is a 3+ dopant with weak magnetic coupling) into 
Lai-ajCa^^MnOs at x = 0.37. 

The left panel in the middle row in Fig. 3 shows the 'flat' 
profile of TT-r in the undoped FM-M state. The central 
panel shows the dopant locations at 77 = 6% and the right 
panel shows the Ur in the presence of dopants. The final 
Ur shows local charge order arising out of 'disordering' the 
FM-M! This has indeed been reported [11] recently. The 
FM order in the ground state is significantly suppressed 
with increasing 7^, and Fig.4.(b) shows the low T metal- 
insulator transition driven by B site doping. Fig.4.(c) 
shows the dramatic suppression in resistivity of this PS 
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Fig. 5: Colour online: Nearest neighbour S^.S^ from MC snap- 
shots in the presence of an applied field. Top row, x = 0.4 with 
2+ dopants at 77 = 0.06, V = 5, = 0. Bottom row, x — 0.5 
with 3+ dopants at ry = 0.06, F = 5, J' = 0.2. The fields are 
(left to right) /i = 0, 0.002, 0.005. The uniform (red) regions 
are FM while the stripelike (green) regions are AF. 



Fig. 4: Colour online: Temperature and field dependence of 
resistivity p(T). (a) Metallisation of the x = 0.25 FM-CO- 
I with non magnetic dopants of valence 2+ and 1/ = 5, the 
curves are for different 77. (b) Ferromagnetic metal to (CO) 
insulator transition at x = 0.40 with increasing concentration 
of non magnetic 2+ dopants with 1/ = 5, (c) Field response 
at X = 0.40 and ry 6%: ho = 0, /ii = 0.002, /i2 0.005. 
(d) Resistivity in the undoped CE-CO-I at x = 0.5 (dotted 
line), 'metallisation' at 77 = 6% with 3+ dopants {V — J' — 
0.2) and field response. 



system, at 77 = 6%, in response to a magnetic field. The 
disorder driven insulating state can be readily 'metallised' 
by a modest field. 

Just as depleting the electron density converts the fer- 
romagnetic metal at x = 0.4 to a FM-M+AF-CO-I phase 
separated state, we can explore the effect of increasing 
the electron density on the CE state at x = 0.5. Since 
n ^ no — r]{3 -\- X — a) it is obvious that dopants like 
Ru or Sn, with valence 4+, will increase the carrier den- 
sity to 71 ^ 0.5 + 0.577. We explicitly checked that this 
leads to suppression of CE order and charge order, en- 
hanced FM correlation, and at some 77c an insulator-metal 
transition [20]. A simpler version of this was explored 
earlier [14]. For B valence 4+ the magnetic character of 
the impurity, or Coulomb interaction Vc, are not relevant. 
However, for B ion valence 3+ the electron density reduces 
to ri ~ 0.5 — 0.57^. In that case the carrier density is not 
in the PS window between CE-CO-I and FM-M but be- 
tween the CE-CO-I and the A-2D phase! The resulting 
state would not be a FM-M unless the magnetic aspect of 
the B dopants affect the background antiferromagnetism. 

We highlight the case where the dopant, e.g^ Cr, has a 
strong AF coupling to the neighbouring Mn. In contrast 
to the reference CE phase with J = 0.1, where a spin has 
two parallel and two antiparallel neighbours, using B ions 



with J' = 0.2 forces all four nearest neighbours (NN) to be 
antiparallel to the B spin. Due to its low Cg occupancy the 
B site acts locally like the 'corner site' of a CE chain, and 
the four next NN spins are also antiparallel to the B spin. 
This '1 + 8' cluster [20] has a large moment but at low r] 
distant clusters are uncorrelated so the magnetic 'recon- 
struction' by itself does not create a global FM state. The 
presence of small NN Coulomb repulsion Vc suppresses the 
Cg charge density on Mn sites which neighbour the B ions. 
To conserve the electron count, the charge that is pushed 
out creates small FM-M regions of locally high density, 
n ~ 0.6, dominated by the double exchange interaction. 
These FM-M droplets connect the otherwise disconnected 
'1 + 8' clusters. The resulting pattern is a complex mix 
of FM-M, FM-CO and AF regions. The ferromagnetism 
in the magnetic B ion doped case emerges from a combi- 
nation of {i) breakup of the CE pattern by the magnetic 
dopants, creating tiny FM clusters, and (ii) their coupling 
via FM-M droplets created by charge pushed out due to 
the Coulomb interaction. 

In addition to converting insulators to metals, and vice- 
versa, B site doping effects are interesting because the fi- 
nal state is necessarily inhomogeneous, with possibly huge 
magnetoresistance (MR). For example if the PS is between 
FM-M and AF-CO-I the randomly located B ions frag- 
ment the PS state and the FM-M domains in such a sit- 
uation are weakly linked to each other at zero field. The 
large 'moments' of the FM-M domains, see Fig. 5, can be 
aligned by a small field leading to enormous increase in 
conductance. This large polarisability and MR are key 
signatures of a PS state. The large changes in resistivity 
driven by an applied field are shown in Fig.4.(c)-(d). Fig. 5 
shows the spatial evolution of the magnetic state for two 
systems doped into PS windows. The dramatic transport 
response correlates with a growth in volume of the FM-M 
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phase. The effect is stronger for the x = 0.4 case, but is 
also visible for the more complex x = 0.5 case with 3+ 
magnetic dopants. 

Since we have attempted to model a complex phe- 
nomenon, and had to make several approximations to 
make progress, let us list out some of the checks [20] we 
have performed to establish the robustness of our results, 
(i) Any real material will have some degree of A site disor- 
der, while we have assumed the reference manganite to be 
'non disordered'. We have checked that our qualitative re- 
sults survive even if we include binary 'A site' disorder of 
magnitude upto 0.1 in our model. In the experiments also 
most B site results are on the La-Pr-Ca family, where the 
cation disorder cr^ is small, (ii) Impurity valence fluctua- 
tion: in reality it is not the valence of the impurity, but its 
electronic level, that is fixed. We had assumed B dopants 
to be of integer valence, but that is true only if the B site 
potential V is large, and the conduction electron density 
rii on the B site is <C 1. We tried out 'softer' B potentials, 
down to y = 1, and that leads to rii ~ 0.25 on the B site. 
For an assumed valence of 3, say, this implies a modest 
change to 2.75, not significantly affecting our argument, 
(iii) Weak localisation (WL) effects: we have described 
'metals' in 2D in the presence of B site disorder. On large 
sizes such a system will show WL effects, and a genuine 
metallic phase will occur only in three dimensions. 

Let us conclude. We have illustrated how one might en- 
gineer phase separation in low cation disorder manganites 
by doping B ions of suitable valence and magnetic char- 
acter. The most promising regimes are close to commen- 
surate filling, X = 0, 0.25, 0.50, etc, which are typically 
ordered insulating states close to a metal. Most of the ef- 
fects can be understood in terms of the valence change of 
Mn, but the B-Mn magnetic coupling is also crucial, par- 
ticularly for dopants on the CE phase. The percolative 
state that emerges is typically highly polarisable and has 
a large low field magnetoresistance. The current paths in 
the phase separated regime are dictated by avoidance of 
B dopant locations. Since presently available techniques 
allow atomic scale manipulation [21] of dopant locations, 
B site doping opens up the prospect of controlling the 
nanoscale current paths in a manganite. 

We acknowledge use of the Beowulf cluster at HRI. 
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